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Abstract 

We consider the possibility of using one of the _D-flat directions in the minimal 
supersymmetric standard model (MSSM) as the inflaton. We show that the flat 
direction consisting of (first generation) left- and right-handed up-squarks as well as 
the up-type Higgs boson may play the role of the inflaton if dominant part of the 
up-quark mass is radiatively generated from supersymmetric loop diagrams. We also 
point out that, if the -R-parity violating Yukawa coupling is of O(10~ 7 ), i?-odd D-flat 
directions may be another possible candidate of the inflaton. Such inflation models 
using D-flat directions in the MSSM are not only testable with collider experiments 
but also advantageous to resolve the problem how the inflaton reheats the universe. 



Inflation pQ is now one of the most important ideas in cosmology. Inflation not only 
solves the serious horizon and flatness problems but also provides a viable scenario of gen- 
erating the origin of cosmic density fluctuations.* 1 In particular, precise measurement of 
the anisotropy of the cosmic microwave background (CMB) by the Wilkinson Microwave 
Anisotropy Probe (WMAP) suggests that the primordial density fluctuations are almost 
purely adiabatic and scale-invariant [3], which are predictions of (some classes of) infla- 
tionary models. In the inflationary models, a scalar field, called "inflaton," is introduced 
to realize the inflationary epoch of the universe. During inflation, potential energy of 
the inflaton gives the dominant part of the energy density. In order to realize the quasi 
de Sitter universe with the potential energy of the inflaton, its kinetic energy should be 
much smaller than the potential energy during inflation. Consequently, we are led to the 
paradigm of slow-roll inflation, where inflation is driven by the potential energy of slowly 
evolving scalar field. 

It is non-trivial to find a scalar field which satisfies the slow-roll condition, and from 
the particle-physics point of view, it is important to find good candidates of the inflaton. 
In particular, it is interesting to ask if the inflaton can be observed at high-energy collider 
experiments in any scenario of inflation. In many classes of models, however, a new scalar 
field is introduced as the inflaton. Usually, such scalar fields do not belong to the standard 
model and have very weak interactions with the standard-model particles, which makes it 
very difficult to find and study the inflaton with collider experiments. Moreover, since the 
interactions between such inflaton and standard-model particles cannot be determined, 
it must be given ad hoc by hand, which obscures the thermal history of the universe. 
Thus the reheating temperature, for instance, cannot surely be estimated. In order to 
construct a testable and economical model of inflation, it is desirable to find a candidate 
of the inflaton in the scalar fields which are in some sense familiar to us. In fact, in the 
standard model, the only scalar field is the Higgs boson. As we will see later, however, 
it is known that the Higgs boson in the standard model cannot be the inflaton since its 
(quartic) coupling is too large to generate cosmic density fluctuations that are consistent 
with the observations. In addition, in the standard model, the potential of the Higgs boson 
is significantly affected by radiative corrections and hence we cannot expect flat enough 
potential which is crucial to cause inflation long enough. Thus, we should conclude that 
it is impossible to find a viable candidate of the inflaton in the standard-model fields. In 
fact, the second point, effects of radiative corrections to the inflaton potential, is in general 
a serious problem in constructing inflation models. 

In order to control the radiative corrections, it is often the case that inflation models are 
considered in super symmetric framework; indeed, in supersymmetric models, quadratic di- 
vergences cancel out between bosonic and fermionic loops and the flatness of the potential 
can be guaranteed. Thus, in considering the inflation, supersymmetry is likely to play very 
important roles, and in this letter, we adopt (low-energy) supersymmetry. If we super- 

Another possibility of generating cosmic density fluctuations may be to adopt the "curvaton" scenario 
P]. Here, we do not consider such a possibility and assume that the cosmic density fluctuations are totally 
generated from the primordial fluctuation of the inflaton. 
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symmetrize the standard model, various scalar particles are introduced as superpartners 
of quarks and leptons. Those scalars may play the role of inflaton. It is worth noting that, 
if this is the case, the reheating processes into the standard-model particles are obvious, 
which makes it possible to study both the inflationary and thermal history of the universe 
only with the low-energy "known" physics. Actually there was an early attempt to build 
inflation models along this line [3]. However, as discussed later, it contained some difficulty 
in producing the right amount of density fluctuations, so the authors of Ref. jl] had to 
introduce additional mini-inflation. Here we would like to pursue another possibility. 

In this letter, we consider the possibility of using scalar fields in the minimal super- 
symmetric standard model (MSSM) as the inflaton (which is denoted as <fi hereafter). In 
particular, we will discuss that the .D-flat direction consisting of first-generation left- and 
right-handed up-type squarks as well as the up-type Higgs boson may play the role of the 
inflaton if the up-quark mass is radiatively generated. We will also see that, if the -R-parity 
violating Yukawa coupling is of O(10 -6 — 10 -7 ), i?-parity violating .D-flat directions may 
also be the inflaton. 

We start with discussing possible scenario of inflation within the MSSM. In the frame- 
work of the slow-roll inflation, the amplitude of the inflaton field varies during and af- 
ter the inflation. The change of the inflaton amplitude is typically of 0(M*), where 
M* ~ 2.4 x 10 18 GeV is the reduced Planck scale.* 2 Thus, inflaton originating from 
the MSSM sector should have an amplitude of 0(M*) during inflation. Consequently, we 
are forced to consider the chaotic inflation [5J within the MSSM assuming that the MSSM 
field corresponding to the inflaton has an amplitude of the order of the Planck scale during 
inflation. Later we will discuss how to realize such a large field value. 

Now, we discuss the observational constraints on the inflaton potential. In the simplest 
scenario, the cosmic density fluctuations are parameterized by the curvature perturbation 
1Z which depends on the inflaton potential V as 



K(k) 



2vr V 



- k=aH ln f 



where H in { is the expansion rate H of the universe during inflation, which is related to the 
potential energy of the inflaton during inflation as H? nf = V/3M%, while V is the derivative 
of the inflaton potential with respect to the inflaton field. Here, notice that this quantity 
is evaluated at the time when the fluctuation exits the horizon during inflation; k and a 
denote the wave-number (for the comoving coordinate) and scale factor, respectively. 

From the measurements of the CMB anisotropies (as well as from other observations), 
we can obtain constraints on TZ(k). The most important constraint is on its normalization. 
Due to the fact that AT/T ~ O(10~ 5 ), the typical size of 1Z is also constrained to be 
O(10" 5 ). More precisely, the WMAP team [3j estimated it as \TZ(k)\ 2 = 2.95 x IQ~ 9 A 

# 2 We consider neither the so-called small-field models (for example, new inflation) nor hybrid models 
within the MSSM. Generally speaking, both require the vacuum to be such that the standard-model gauge 
symmetries are spontaneously broken, which is hardly justified. Here and hereafter we concentrate on the 
so-called large-field models (i.e., chaotic inflation). 
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with A = 0.9 ± 0.1 at k = 0.05 Mpc -1 , assuming power-law ACDM model. If the inflaton 
potential has the parabolic form V = |M|0 2 , the inflaton mass is required to be ~ 
O(10 13 GeV). Obviously, such a heavy scalar field does not exist in the MSSM. 
We consider the next possibility where the inflaton potential is quartic: 

V = \\<P\ (2) 

where A is a dimensionless coupling constant. Then, the normalization of the primordial 
density fluctuation requires A ~ O(10~ 13 ) as follows. The value of the inflaton field is 
related to the e-folding number N as <fi ~ \/8NM*. (See Eq. (jl4j) below). Thus the 
curvature perturbation at the horizon exit is evaluated as 1Z ~ 0.3-\/AA^ 3 / 2 . Equating this 
with the WMAP result, we obtain A ~ 10 -13 , where iV ~ 60 is used. Before studying the 
MSSM cases, we should comment on the non-supersymmetric case; from this constraint, 
it is obvious that the Higgs boson cannot play the role of the inflaton since, if the quartic 
coupling of the Higgs boson is as small as O(10~ 13 ), Higgs mass becomes O(10~ 4 GeV) 
which is unacceptably smaller than the present experimental bound. 

In the supersymmetric case, such a small coupling for quartic interaction cannot be 
realized if the potential is lifted by the gauge D-term interactions, since, if so, the coupling 
constant A becomes of the order 0(g 2 ) where g is the gauge coupling constant in the 
standard model. Therefore we focus our attention on the D-flat directions. For D-flat 
directions, we have to be more careful since behaviors of the potential depend on which 
flat direction we consider. In the MSSM, Yukawa interactions exist in the superpotential to 
generate the fermion masses. Such Yukawa interactions lift some of the D-flat directions. 
In addition, we can also find several D-flat directions which are not affected by the Yukawa 
interactions associated with the fermion masses; without .R-parity violation, such D-flat 
directions are only lifted by the effects of supersymmetry breaking.* 3 (See Ref. jH] for the 
details.) 

We first consider the D-flat direction lifted by the .R-parity conserving Yukawa interac- 
tions. D-flat directions are parameterized by gauge invariant monomial of the superfields. 
We denote the MSSM superfields as Q(3,2,|), Z7(3*, 1, — |), D(3*, 1, §), L(l,2,-|), 
E(l, 1, 1), H u (l,2,^), and H d (l, 2, — |), where we show the quantum numbers for the 
SU(3)c x SU{2)l x U{1)y gauge group in the parentheses. Then, the relevant part of the 
MSSM superpotential is given by 

W = [Yu^QiUjHv + \Y D ]ijQiDjH d + [Ye^L^H,, (3) 

where i and j are generation indices. If we consider the D-flat directions represented 
as QiUjH u , QiDjH d , and LiEjH d , those D-flat directions acquire quartic potential due 
to the Yukawa interactions, as given in Eq. (J2J). In order to relate the Yukawa coupling 
constants to A in Eq. (J2J), we express those D-flat directions with a complex scalar field 

# 3 Here, we assume that coefficients of non-renormalizable terms are suppressed enough to be neglected. 
This may be explained by the f?-symmetry, assigning i?-charge | to each MSSM chiral superfields. 
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$, for example, 

Q ' = 7S&' = > ^7s(*)' (4) 

Self quartic coupling constants of those flat directions, denoted as Xq^Hui ^Q^Ha, and 
^LiEjHd, respectively, are then given by 

XQ t u 3 H u = - \ [Yu]ij\ 2 , X QiDjHd = - \ [YD]ij\ 2 , X LiE] H d = - \ [YE}ij\ 2 , (5) 

where we defined </> = V^Re $. Thus, if one of these coupling constants is of 0(1CT 13 ), we 
may have inflaton candidates within the MSSM particles. 

If the fermion masses totally originate from the superpotential given in Eq. Yukawa 
coupling constants are evaluated as 

[Yu}u = — r^, [Yd\u = [Ye]u = (6) 

vsmp v cos [3 v cos p 

where Yjj, Yd, and Ye are Yukawa matrices in the diagonalized basis, while m Ui , m di , 
and m ei are fermion masses for corresponding Yukawa couplings. In addition, v = 



(H%) + (H%) ~ 174 GeV while tan/3 = (H u )/(H d ). Yukawa coupling constants for 
the second and third generation quarks and leptons are much larger than ~ 10~ 6 , and 
hence we consider the possibility of using the first generation squarks and/or sleptons as 
the inflaton. (Hereafter, we consider only the first generation squarks and sleptons, and 
drop the generation indices for simplicity unless otherwise mentioned.) If we estimate the 
Yukawa coupling constants for the up and down-quarks as well as electron using Eq. (|Hjl. 
we obtain 

8.6 x 1(T 6 / m u \ 2.9 x 1(T 5 / m d \ 2.9 x 10" 6 

Vu ~ sin/5 X U.5 MeV J ' Vd ~ cos (3 X \5MeV ) ' Ve ~ cos (3 ' ( ' 

In the MSSM, it is often the case that, in order to evade the Higgs-mass constraint, 
relatively large value of tan (3 is required [7j. Then, y d and y e are likely to be larger 
than 10~ 6 . In fact, if we adopt the up-quark mass m u = 1.5 — 4.5 MeV [8J,* 4 we obtain 
Vu > 8.6 x 10~ 6 , leading to \quh u > 2.5 x 10~ n , and hence even the up-squark cannot 
play the role of the inflaton.* 5 One might well give up the simple A</> 4 model and add 
another mini-inflation to reconcile the predicted magnitude of density fluctuations with 
the observed one jl]. However, here we would like to stick to the A</> 4 model, since it 
predicts the density fluctuations with very small uncertainty as shown later. 

# 4 Here, we do not consider renormalization group running which suppresses the fermion masses at 
higher energy scale. Even if the renormalization group effects are taken into account, our discussions are 
qualitatively unchanged. 

# 5 The quark mass ratios are constrained by imposing a limit on next-to-leading order corrections in 
the chiral perturbation theory. If the corrections become sizable, the up quark mass can be much smaller 
than the lower limit ~ 1.5 MeV Then, y u can be as small as ~ 0(1CP 7 ). In this case, we do not have 
to consider the radiatively induced up-quark mass as we will discuss in the following. 
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Figure 1: Loop diagram contributing to the up-quark mass in the MSSM. Black dots 
represent the insertion of the off-diagonal elements of the squark mass matrices while the 
vertex with the open circle is from the insertion of the trilinear coupling A t . 



So far, we have assumed that the Yukawa interactions in the superpotential are the 
only sources of the fermion masses. In the MSSM, however, radiative corrections due 
to the supersymmetric loop also affect the fermion masses. If such an effect is large 
enough to explain the dominant part of the fermion mass, Yukawa coupling constant in 
the superpotential can be smaller than the naive expectation. Indeed, for the up-quark 
mass, contribution from the squark-gluino loop diagram may become sizable with non- 
vanishing off-diagonal elements of the left- and right-handed squark mass matrices, as 
given in Fig. ^ Relevant part of the supersymmetry breaking terms contributing to this 
diagram is given by 



C = Aml L - tL u L t* L + Aml RiR u R t* R + A t t L t R H u + ^m^GG + h.c, 



where Ul,r an d t^R are left- and right-handed up- and top-squarks, respectively, while 
G denotes the gluino. Hereafter, we approximate that the masses of the squarks are all 
degenerate for simplicity, and we define 



xs % &* s !<£±, a ,^, (9) 



where rrig is the squark mass, and y t is the top Yukawa coupling constant. Then, the loop 
contribution to the up-quark mass is given by 

™f> = ^dm^\5>W - fa | +_^ + 2 + telnl » , (10) 

where g 3 is the SU(3)c gauge coupling constant. 

Size of the off-diagonal elements of the squark mass matrices are constrained from 
flavor changing processes. In particular, for one finds ^ < 4.6 x 10~ 2 (9.8 x 1(T 2 , 
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2.3 x 1CT 1 ) for x = 0.3 (1.0 and 4.0) and m- q = 500 GeV [Ej. #6 Constraint on 5\f> is, 
on the other hand, not available since we do not have precise measurement of the flavor 
changing decay of the top quark. Thus, the combination 5[^5[^ can be of O(0. 1) without 
conflicting experimental constraints.* 7 Then, for example for the case where mg ~ rriq, 
we obtain 

m ^°° P) - ^sWi?^?^) * 3.6 MeV x y t a t sm(3 I J f^j • (H) 

Thus, the up-quark mass of the size of a few MeV can be generated from the loop effect 
with reasonable value of a t ~ 0(1). Then, the up- Yukawa coupling constant in the 
superpotential can be of O(10 -6 — 10~ 7 ). Notice that the loop-induced up-quark mass 
is not suppressed even if the masses of the superparticles become much larger than the 
electroweak scale (as far as all the supersymmetry breaking parameters are of the same 
order). Thus, if one wishes to push up m^° op ^ without affecting flavor- changing processes, 
one possibility is to assume (relatively) large masses for the superparticles. We can also 
consider radiative correction to the down quark mass. In this case, A t should be replace 
by Vbl^H in the calculation, where y& and \lh are bottom-quark Yukawa coupling and 
super symmetric Higgs mass, respectively, and the generation mixings are provided by 
the down-type squark mass matrices. Thus, radiative correction to the down-quark mass 
may become sizable when tan/? is large. Constraints on the off-diagonal elements of the 
down-type scalar quark mass matrices are, however, more stringent than those on the 
up-type scalar quark mass matrices. Thus, it is more difficult to radiatively generate the 
down-quark mass. 

The above arguments open a window to use the D-flat direction QUH U as the inflaton. 
Thus, hereafter, we concentrate on the evolution of the universe for the case where D-flat 
direction QUH U plays the role of the inflaton, assuming that the dominant part of the 
up-quark mass is from the radiative correction. 

The potential of QUH U is now given by 

V&) = \rriy + (12) 

where ~ O(100 GeV) is the soft supersymmetry breaking mass. These two terms 
become comparable at (fi m = \j2/\ m§. Since the mass term can be neglected for large 
value of ^> (f) m , its dynamics is just the same as the well-known A0 4 model. First let 

# 6 Constraint on here is from the mass difference of the neutral Bd mesons. Such a constraint is 
obtained from the basis where the down-type squarks are diagonalizcd. Two basis are related by the 
Kobayashi-Maskawa matrix, but we neglect the difference of the two basis because the change of the 
constraint is small enough to be neglected. 

# 7 If the supersymmetry breaking parameters have large complex phases, too large neutron electric dipolc 
moment may be induced in this model. Here, we consider the case where the phases in the supersymmetry 
breaking parameters are small enough. We also assume that the mixing parameters for the first and second 
generation squarks are suppressed in order to evade the constraints from the K°-K° mixing. 
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us briefly review the A0 4 model and discuss the observational constraint on that. In this 
model, the slow roll parameters are given by 




Ml , I '" Ml , . 

= V "" M *V = V' (13) 

where we substituted the quartic potential Eq. (J2J). When the inflaton becomes equal 
to cn d = 2\/3M*, the slow-roll condition breaks down, leading to the end of inflationary 
epoch. The e-folding number N is given by 

/■tend r<t> 1 

where we used the slow-roll approximation in the second equality. 

In order to impose observational constraints on inflation models, it is necessary to 
evaluate several parameters that characterize the density fluctuations generated during 
inflation. They are scalar spectral index n s , its running dn s /d\nk, and tensor-to-scalar 
ratio r, apart from the normalization that determines the value of A. These are related to 
the e-folding number as 

- 1 3 dn s 3 _ 16 

N + l' d\nk~ r ~N+l 1 J 

Thus we must precisely evaluate N in order to make a comparison between predictions 
and observational results. Since the value of N necessary to solve the horizon and flatness 
problems depends on the thermal history of the universe, we must specify the reheating 
processes. Fortunately, however, it is possible to determine the e-folding number as N ~ 64 
with small uncertainty in the case of \(j) 4 model. The reason is that the energy density 
of inflaton oscillation decreases just like radiation after the end of inflation, which applies 
only to the quartic model. The explicit expression for N is given by ^1] 

^ = 62-ln^)-ln(^ 7TT ^)+ln(^ 7J ], (16) 

where V represents the potential energy of the inflaton, and the subscripts "0," u k," 
and "end" are for variables at present, at the time when the fluctuation with the wave- 
number k exits the horizon, and at the end of the inflation, respectively. The right hand 
side of Eq. (jT6|) gives N ~ 64 for k = aoH , irrespective of the details of the reheating. 
Here we have assumed that the reheating processes complete before the amplitude of the 
inflaton becomes smaller than <fi m , which will be justified below. Substituting A^ ~ 64 into 
Eqs. (115(1 . we obtain n s ~ 0.95 and r ~ 0.24 with negligible running of n s . 

Constraints on the inflationary models driven by a single slow-rolling scalar field from 
the recent observations including the WMAP have been studied extensively ^21 El E] • 
Although Ref. ^2] claimed that A0 4 model lies in the region marginally excluded by the 
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WMAP data in combination with smaller scale CMB and large scale structure survey 
data, more detailed analyses [T3J Ej showed that the model cannot be excluded by the 
WMAP data alone for N^AO. In addition, the recent systematic study using both the 
Sloan Digital Sky Survey and WMAP demonstrated that A0 4 model is still allowed [TK] . 
Thus there is no reason to disregard A0 4 model at present, and one could well contend 
that the model is of much interest since it is on the edge. 

The next discussion concerns the decay processes of the inflaton. After inflation ends, 
the inflaton oscillates around its origin until the decay completes. In the usual chaotic in- 
flation models, the reheating process proceeds through nonperturbative particle creation 
(preheating) [TH]. During the oscillation, the particles coupled to the inflaton are pro- 
duced when the so-called adiabatic condition is violated, i.e., lu/uj 2 ^!, where uo is the 
effective frequency of the produced particles. In our model, the inflaton field is complex, 
and its nontrivial trajectory in the potential may lessen the efficiency of the preheating 
process ^7]. The possible source of the nontrivial orbit of the inflaton in our case is the 
A-term such as 

V A = a u y u mq$ 3 + h.c, (17) 

where a u is a constant of 0(1). The effect of this A-term is, however, so small that the 
inflaton $ exhibits almost straight-line motion on the complex plane. Therefore, when 
the inflaton comes closest to the origin, its amplitude is much smaller than the critical 
value, below which the adiabaticity condition is violated. Hence the preheating should 
proceed in the same way as a real scalar field. In fact we have confirmed numerically that 
the instability band almost coincides with that in the case of a real scalar field, even in 
the presence of the A-term. Also, since the preheating occurs very efficiently and ceases 
within several oscillations as shown below, cumulative disturbance of the homogeneous 
motion caused by the A-term can be neglected. 

We would like to focus on the four-point scalar interaction with the stop among many 
interactions the inflaton feels. When the inflaton first reaches = 0, the particle produc- 
tion occurs and the stops are generated, typically with the momentum k Tes ~ y\ A 1 / 4 ^ 
and the occupation number nu ~ 0(1) for k ~ k res , where 0o is the amplitude of the 
oscillation. Since the stop obtains an effective mass m 4 ~ ~ yt\4>\ through the four-point in- 
teraction, generated stops are fattened as the amplitude of inflaton increases, and become 
as heavy as ~ ytM* around the endpoint of the oscillation of the inflaton. The stop can 
then decay into two fermions, wino and (left-handed) bottom quark, for example, because 
the time scale of the decay is much shorter than that of the oscillation of the inflaton. 
Similar process occurs when reaches = again. Each time passes its origin, stops 
are generated, and they decay when reaches around its maximal value. This type of 
preheating is known as "instant preheating" ^Hj- Applying the result of Ref. ^Hj to our 
case, the amount of the energy dissipated by each decay of the stops is comparable to that 
originally stored in the inflaton. Thus, the decay processes of the inflaton proceed very 
efficiently and complete within several oscillations. 

After that, the decay products are quickly thermalized through decays, scatterings, 
and annihilations by gauge interactions, which proceed very efficiently. The reheating 
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temperature Trh is expected to be very high: Trh ~ O.IA^M* ~ 10 14 GeV. Such high 
reheating temperature might lead to the overproduction of dangerous relics like graviti- 
nos [T9] . In order to avoid this problem, we assume that one of the followings is realized. 
One solution is to assume large late-time entropy production from, for example, thermal 
inflation [20J. Another is to have a relatively heavy gravitino mass, m 3 / 2 — 10 — lOOTeV, 
so that the gravitinos can decay well before the big bang nucleosynthesis (BBN) epoch, 
Xg3N ~ lMeV. However, this does not remedy the situation if the lightest supersymmet- 
ric particle (LSP) is the standard bino-like LSP, since those produced from the decay of 
gravitinos would overdose the universe. This difficulty can be evaded by the introduction 
of a supersymmetric partner with a mass much lighter than lOOGeV. One possibility may 
be the axino, superpartner of the axion |2T] . 

So far we have not mentioned how the flat direction can take the value of 0(M*) or 
larger avoiding the Hubble-induced mass term that prevents the flat direction from slow- 
rolling (so called 77-problem). In the minimal supergravity model, the scalar potential 
includes an exponential factor which essentially precludes any scalar amplitudes larger 
than M*, and scalar masses of the order of H are generated. An idea to circumvent these 
obstacles to construct a successful chaotic inflation model is to introduce the Heisenberg 
symmetry [22] under which the inflaton and a chiral field z transform as follows, 

5z = e*0, 8<j) = e, (18) 

where e is a complex parameter. We can construct an invariant combination y from z and 

as 

y = z + z* -0*0. (19) 

Imposing the Heisenberg symmetry in the Kahler potential, it is written only with y, i.e., 
K — f(y). It is easy to see that y and should be regarded as independent variables 
since the kinetic terms are diagonalized for these variables [2B|. It was shown that this 
symmetry protects the flatness of the inflaton potential from both the exponential growth 
and Hubble-induced mass term. It should be emphasized that the introduction of a new 
degree of freedom y is the price for keeping the potential flat. Central to this issue is the 
problem how to stabilize y. The dynamics of y with a specific form of f(y) = |lny + y 2 
was discussed in Ref. [21], and it was found that the value of y is fixed during the inflation. 
In addition, another mechanism of fixing y using the radiative corrections was proposed 
in Ref. [2Hj i n the case of the no-scale supergravity model (2H| where the Kahler potential 
takes a special form as K = —3 Iny. In both cases, the potential of the inflaton is same as 
that in the global supersymmetry case as long as y is fixed. In this letter, we just assume 
y is somehow fixed and remains constant for successful chaotic inflation.* 8 

# 8 Notice that, even in this framework, soft supersymmetry breaking terms required for our mechanism 
can be generated by introducing a new supersymmetry breaking field x (with a slight modification of the 
Kahler potential). For example, let us consider the Kahler potential of the form K = f(y) + aij \x\ 2 (cf>* 4>j + 
h.c.) with a,ij being constants, 4>i MSSM chiral multiplets, and y = z + z* + (f)*4>i + \x\ 2 . This Kahler 
potential does not have the Heisenberg symmetry if ^ 0. Previous arguments, however, still hold since 
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Finally, we comment on the case with the .R-parity violation. In the MSSM, there are 
D-flat directions which are not lifted by the Yukawa interactions given in Eq. Q. Such flat 
directions are parameterized by the monomial of the superfields with odd .R-parity, and can 
be lifted if the -R-parity is broken. Choosing relevant .D-flat direction, .R-parity breaking 
Yukawa coupling can be as large as ~ O(10~ 6 — 1CT 7 ) without conflicting experimental 
bounds if the baryogenesis takes place after the sphaleron interaction becomes ineffective 
|26j . Thus, if we adopt -R-parity violation of this size, such an .R-parity violating .D-flat 
direction can be another candidate of the inflaton. The decay of .R-parity violating .D-flat 
directions occurs in the same way as the previous case via the "instant preheating." Since 
they do not necessarily have interactions with stops, the four-point scalar interaction with 
the largest coupling constant is expected to come from the .D-term. In particular, there 
are four-point scalar interactions ~ g 2 x 2 4> 2 , where x represents the field orthogonal to the 
flat direction. 

In summary, we have investigated whether the inflation can be embedded in the MSSM 
sector, and found that the D-flat direction consisting of the first generation left- and right- 
handed up squarks and the up-type Higgs boson may be the inflaton if the up-quark mass 
predominantly comes from the one-loop threshold correction to the up- Yukawa coupling 
constant. The dynamics of the inflaton is almost the same as the A0 4 model, which has 
attracted much attention recently. Since the inflaton in our model consists of the MSSM 
particles, it is not only minimal but also testable at high-energy collider experiments. In 
particular, if the QUH U flat direction plays the role of the inflaton, mixings of the up- 
and top-squarks should be large. This is an interesting check point of our model and 
can be tested by collider experiments as well as precise measurements of flavor- changing 
processes. 
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